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Abstract 

A search for TeV - PeV muon neutrinos from unresolved sources was performed on AMANDA-II 
data collected between 2000 and 2003 with an equivalent livetime of 807 days. This diffuse analysis 
sought to find an extraterrestrial neutrino flux from sources with non-thermal components. The 
signal is expected to have a harder spectrum than the atmospheric muon and neutrino backgrounds. 
Since no excess of events was seen in the data over the expected background, an upper limit of 
£'^$go%c.L. < 7.4 X 10~^ GeV cm~^ s~^ sr~^ is placed on the diffuse flux of muon neutrinos with 
a $ oc spectrum in the energy range 16 TeV to 2.5 PeV. This is currently the most sensitive 
$ oc diffuse astrophysical neutrino limit. We also set upper limits for astrophysical and prompt 
neutrino models, all of which have spectra different than $ oc E~^. 

PACS numbers: 95.55.Vj, 95.75.-z, 95.85.Ry 



* Corresponding author: hodges@icecube.wisc.edu (J. Hodges) 

ton leave of absence from Universitat Erlangen-Niirnberg, Physikalisches Institut, D-91058, Erlangen, 
Germany 

ton leave of absence from Universita di Bari, Dipartimcnto di Fisica, 1-70126, Bari, Italy 

^affiliated with Dept. of Chemistry and Biomedical Sciences, Kalmar University, S-39182 Kalmar, Sweden 

^NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 

"affiliated with IFIC (CSIC-Universitat de Valencia), A. C. 22085, 46071 Valencia, Spain 



5 



I. INTRODUCTION 



High energy photons have been used to paint a picture of the non-thermal Universe, but 
a more complete image of the hot and dense regions of space can potentially be obtained 
by combining astrophysical neutrino and gamma ray data. Neutrinos can provide valuable 
information because they are undeflected by magnetic fields and hence their paths point back 
to the particle's source. Unlike photons, neutrinos are only rarely absorbed when traveling 
through matter. However, their low interaction cross section also makes their detection 
more challenging. The observation of astrophysical neutrinos would confirm the predictions 

, |9| that neutrinos are produced in hadronic interactions in cosmic 
accelerators, such as active galactic nuclei or gamma-ray bursts. 

Instead of searching for neutrinos from either a specific time or location in the sky, this 
analysis searches for extraterrestrial neutrinos from unresolved sources. If the neutrino flux 
from an individual source is too small to be detected by point source search techniques, 
it is nevertheless possible that many sources, isotropically distributed throughout the Uni- 
verse, could combine to make a detectable signal. An excess of events over the expected 
atmospheric neutrino background would be indicative of an extraterrestrial neutrino flux. 

In this paper, we report on a search for a diffuse flux of astrophysical muon neutrinos 
performed with data collected by the AMANDA-II neutrino telescope from 2000 - 2003. To 
perform the search, a 5.2 sr sky region (slightly less than 27r sr) was monitored over a four 
year period, for a total of 807 days of livetime. Before describing specifics of the analysis, 
the existing diffuse neutrino models and limits and how we aim to detect neutrinos are 
described in Sections Hi] and IIIII In Section IIV[ typical backgrounds to the extraterrestrial 
signal are discussed, as well as how events are simulated in the detector. We also explain how 
an atmospheric neutrino sample was obtained. An extensive systematic uncertainty study 
is described in Section |Vl The relationship between up and downgoing events is explored 
in Section IVII Finally, the results of the analysis are presented in Section I VI II Since no 
excess of high energy events was seen above the predicted atmospheric neutrino background, 
we set limits on the flux of extraterrestrial muon neutrinos with a generic $ oc energy 
spectrum as well as with a number of different model spectra discussed in Section IVIIBI 
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II. ASTROPHYSICAL FLUXES AND LIMITS 



The analysis presented in this paper assumes a $ oc spectrum resulting from shock 
acceleration processes. Although other spectra were tested, the $ oc spectral shape is 
considered a benchmark to characterize acceleration in many sources. 

The Waxman-Bahcall upper bound follows an $ oc E ^ spectrum and reaching 

below the sensitivity of this bound has traditionally been a goal of neutrino experiments. 
Nellen, Mannheim and Biermann and Becker, Biermann and Rhode js] have suggested 
$ oc E~^ neutrino spectra with higher normalizations than the Waxman-Bahcall bound. 
The other astrophysical neutrino models tested here (Mannheim, Protheroe and Rachen , 
Stecker, Done, Salamon and Sommers {3, Q] and Loeb and Waxman 9]) predict different 
spectral shapes and are specific predictions of neutrino fluxes from classes of objects such as 
active galactic nuclei (AGN) and starburst galaxies. The models have been derived based on 
a variety of astronomical results, including the observed extragalactic cosmic ray flux and 
x-ray and radio measurements. 

A precursor to this muon neutrino analysis was conducted with data collected in 1997 by 
the AMANDA-BIO detector (In 1997, the AMANDA detector consisted of 10 sensor 
strings, a subset of the 19 strings in the final AMANDA-II configuration.) Other AMANDA 
analyses have focused on the search for a diffuse flux of neutrinos using particle showers or 
cascades ll|. Cascades are caused by Uf, and u^- charged current interactions and all- flavor 



neutral current interactions in the ice near the detector. Even though no extraterrestrial 
signal has been detected, models can be excluded by setting upper limits. 



The Frejus 



121, MACRO 



13| and Baikal |l^) experiments have set upper limits on the 



flux of astrophysical neutrinos in the same energy region as this analysis (TeV - PeV). Pub- 
lished upper limits from these experiments assuming a $ oc E^^ spectrum are summarized 
along with the results of this analysis in Section IVII A[ Depending on the detector and 
the specific analysis, the reported upper limit constrains either the muon neutrino flux or 
the all-flavor neutrino flux. Upper limits obtained from all-flavor analyses are not directly 
comparable to upper limits. However, for a wide range of neutrino production mode 



and oscillation parameters, the flavor flux ratio at Earth can be approximated as 1:1:1 15|. 



In that case, either a single-flavor limit can be multiplied by three and compared to an all- 
flavor result, or an all- flavor limit can be divided by three and compared to a single-flavor 
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result. 

The Baikal experiment has placed limits on models with spectra other than $ oc E^^ jl^ , 
which are compared to the results from this analysis in Section IVIIBI Here, nine different 
spectral shapes are tested, including the search for prompt neutrinos from the decay of 
charmed particles. Since this analysis is optimized on energy- dependent parameters, the 
optimization was performed individually for each energy spectrum. 



III. NEUTRINO DETECTION IN AMANDA 



Although chargeless particles like neutrinos are not directly observable, the by-products 
of their interactions with polar ice or rock near the detector can be observed. In particular, 
two types of neutrino-induced events can be distinguished in AMANDA. All neutrino flavors 
can cause hadronic or electromagnetic showers in the ice and these appear as a momentary 
point-like source of Cherenkov light. Alternatively, long track-like events are the signature 
of neutrino-induced muons traveling through the detector. A cone of Cherenkov light is 
emitted by these muons as they travel faster than the speed of light in ice. The present 
analysis focuses exclusively on the muon track channel for identifying neutrino events. Tau 
neutrinos can undergo charged current interactions and contribute to the upgoing /i and z/^ 
fluxes via r — > fiUri^ti decay. Although interactions and r decay may contribute between 
10% to 16% to the signal flux [16], this contribution is ignored in this analysis. 

Nineteen vertical strings hold the optical modules (OMs) for recording the timing and 
position of detected photons, which is needed for reconstructing the path of the muon \v}\. 
The angular distribution between the neutrino direction and the reconstructed muon track 
has a median of 2° when the highest quality events are used. The 677 OMs each consist of 
a photomultiplier tube (PMT) enclosed in a pressure-resistant glass sphere. The OMs are 
deployed to depths between 1500 and 2000 meters. An event is recorded when at least 24 
OMs report seeing light within a 2.5 //s window. AMANDA has been operating in the final 
configuration with 19 strings (AMANDA-II) since 2000 0. 
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IV. SEARCH METHODS 



This analysis uses the Earth as a filter to search for upgoing astrophysical neutrino- 
induced events. The background for the analysis consists of atmospheric muons and neutri- 
nos created when cosmic rays interact with Earth's atmosphere. The majority of the events 
registered in the detector are atmospheric muons traveling downward through the ice. 

Conventional atmospheric neutrinos arise from the decay of pions and kaons created 
in cosmic ray interactions with the atmosphere. Atmospheric neutrinos are able to travel 
undisturbed through the Earth. They can be separated from atmospheric muons by their 
direction, namely by demanding that the reconstructed track is upgoing. The conventional 
atmospheric neutrino flux asymptotically approaches a $ ~ E^^'^ spectrum in the multi- 
TeV range. Prompt neutrinos are the counterpart of the conventional atmospheric neutrino 
flux and will be discussed in Section IVIII 

In the initial sample of 5.2 x 10^ events, many downgoing events were misreconstructed 
as upgoing tracks. Misreconstruction happens because photons scatter in the ice, causing 
directional and timing information to be lost. Hence, the selected upgoing event sample 
not only contains truly upgoing neutrinos, but a certain fraction of downgoing atmospheric 
muons. 

An energy-correlated observable was used to separate neutrino-induced events since the 
predicted astrophysical neutrino flux has a much harder energy spectrum ($ oc E"^) than 
the conventional atmospheric neutrinos from pions and kaons. Any excess of events at high 
energy over the expected atmospheric neutrino background indicates the presence of a signal. 

The search method can be summarized by the following three selection steps: 

1. Use the zenith angle from the reconstructed track to reject obviously downgoing events. 

2. Select events that have observables more consistent with typical long upgoing tracks. 
This separates truly upgoing events from misreconstructed downgoing events. 

3. Use an energy-related observable (number of OMs triggered) to separate upgoing at- 
mospheric neutrinos from upgoing astrophysical neutrinos. 

This analysis relied on simulated data sets of background and signal events. Sixty-three 
days of downgoing atmospheric muons were simulated with CORSIKA 



version 6.030 and 

the QGSJETOl hadronic interaction model. The events were simulated with a $ oc E^^ '' 



primary energy spectrum. These downgoing events are so frequent (~80 Hz at trigger 
level) that two atmospheric muon events produced by unrelated primaries often occur in the 
detector during the same detector trigger window of 2.5 fis. Timing patterns of the light 
from the two tracks may be such that the reconstruction results in a single upgoing track. 
These coincident muon events may be caused by two muons which are each individually 
incapable of triggering the detector with at least 24 OM hits. However, events which only 
hit a few OMs can now trigger the detector when in coincidence with another event. This 
means that a simple trigger rate calculation of 80 Hz x 80 Hz x 2.5 fis is not possible since 
all combinations of events with a total of at least 24 hits can trigger the detector. These 
coincident muon events were simulated for 826 days of livetime and have a frequency of 
about ~2-3 Hz at trigger level. 

Muon neutrinos with a $ oc spectrum were simulated with nusim [lol and reweighted 



to atmospheric neutrino flux predictions [20|, 121|, 122|, 123|, 



24| . as well as to an astrophysical 



muon neutrino flux of E^<l> = 1x10^^ GeV cm~^s~^sr~^. The normalization of the test signal 
spectrum, which is irrelevant when setting a limit, was taken to be approximately equal the 
previous upper limit from the AMANDA-BIO diffuse analysis jlo| . 



A. Filtering the Data Set 

The 2000 - 2003 analysis covers 807 days of detector livetime between February 2000 and 
November 2003. Because of summer maintenance operations, no data were used from the 
polar summer seasons. In the first stage of the analysis, reconstruction software was used to 
make an initial hypothesis on the track direction of every event based on the timing pattern 
of the detected light [l7|. 

Figure [T] shows the zenith angle of the reconstructed tracks for all events at the beginning 
of the analysis (level 0). Vertically downgoing tracks have a reconstructed zenith angle of 0° 
(cos(^)=l). The data set was reduced to 8 x 10^ events by removing all events reconstructed 
with zenith angles less than 80° (cos(^)=0.17). The remaining data set mainly consists of 
misreconstructed downgoing muons and events near the horizon. 

The reduction of the data by three orders of magnitude with the simple zenith requirement 
made it feasible to perform more CPU-intensive track reconstructions on the remaining 
events. Track parameters were adjusted to maximize the log likelihood, given the observed 
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FIG. 1: The cosine of the reconstructed zenith angle is shown for every event at the beginning 
of the analysis (level 0). The experimental data is dominated by downgoing atmospheric muons. 
Events reconstructed as upgoing appear on the left side of the plot and downgoing events appear 
on the right. 

light pattern. Many of the Cherenkov photons scatter multiple times as they travel through 
the ice and this changes their direction and delays the times at which they are likely to be 
detected. An iterative technique was performed in which each event was reconstructed 32 



times 171], each time with a different seed. Each iteration shifts the zenith and azimuth of 
the track and moves the track to pass through the center of gravity of the hits. The best 
track found by the iterative search was used throughout the later stages of the analysis. 

In order to prevent any inadvertent tuning of the event selection criteria that would bias 
the result, a blindness procedure was followed which required that further event selections 
were developed only on simulation and low energy data, where the signal is negligible com- 
pared to the background. The number of OMs triggered (from now on indicated by N^h, 
or number of channels hit) is the energy-correlated observable used to separate atmospheric 
neutrinos from astrophysical ones (Figure [2]). Only low energy data events (low A^ch val- 
ues) were compared to simulation. High energy data events (high A'ch values) were only 
revealed once the final event selection was established. Energy and N^h are correlated since 
high energy events release more energy in the detector causing more hits than low energy 
ones. However, the correlation is not perfect since high energy events occurring far from the 
detector may trigger only a few OMs. 

Event selection was based on observables associated with the reconstructed tracks llTi and 
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FIG. 2: The number of OMs hit during an event (A'ch) was used as an energy-correlated observable. 
Each line on this A^^^h distribution represents events with approximately the same simulated energy. 
High energy events may not be contained within the detector and hence can trigger a wide N^h 
span. 

is described in more detail in Appendix A. In order to separate misreconstructed downgoing 
events and coincident muons from the atmospheric and astrophysical neutrinos, events were 
required to have observables consistent with long tracks and many photons with arrival 
times close to those predicted for un-scattered propagation. The number of photons arriving 
between —15 and +75 ns of their expected un-scattered photon arrival time is referred to 
as the number of direct hits (N^n)- The direct length (Ljir) is the maximum separation 
of two direct hits along the reconstructed track. The smoothness (S) is a measurement 
of how uniformly all hits are distributed along the track and it varies between —1.0 and 
1.0. Positive values of the smoothness indicate more hits at the beginning of a track and 
negative values indicate more hits occur toward the end. Evenly distributed hits will have 
smoothness values near 0. The median resolution (MR) is calculated from a paraboloid fit to 
the likelihood minimum for the track 13]. This method analyzes the angular resolution on 
an event-by-event basis. Lastly, high quality events have higher values of the logarithm of the 
up-to-down likelihood ratio, AL = (— log£down) — (— log/^up)- The likelihoods C^p and £down 
are the product of the values of the probability density function for the observed photon 
arrival times, for the best upgoing and zenith- weighted downgoing track reconstruction jl^ , 
respectively. A more strict requirement for the likelihood ratio was applied to vertical events 
than for events near the horizon. Horizontal events tend to have smaller likelihood ratios 
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since the zenith angle difference between the best upgoing and zenith-weighted downgoing 
track hypothesis is often small. 

The event selection requirements were successively tightened, based on the reconstructed 
track parameters, establishing five quality levels. At Level 5, used for the final stages 
of the analysis, the event sample is expected to contain only truly upgoing events. The 
zenith angle distribution for the events at each quality level is shown in Figure [31 Although 
the entire upgoing zenith angle region is being studied, the event selection requirements 
preferentially retain vertically upgoing events. Horizontal and vertical events must pass the 
same requirements for track length and number of direct hits, however this is more difficult 
for horizontal events since the detector is not as wide as it is tall. 



B. Separating Atmospheric Neutrinos from Astrophysical Neutrinos 

Figure H] shows the Nch distribution for events at Level 5. The optimal place for the 
energy-correlated event observable requirement was established with the simulation by min- 



imizing the expected Model Rejection Factor (MRF) 26|]. The Feldman-Cousins method 



was used to calculate the median upper limit [27|]. The MRF is the median upper limit 
divided by the number of predicted signal events for the z/^ signal being tested. The MRF 
was calculated for every possible A'^^h value and was at its minimum at A'ch > 100. Hence, 
the optimal separation of astrophysical and atmospheric neutrinos is achieved with this A^ch 
requirement. 

The final event sample was composed of events which pass all event selection requirements 
(Level 5) and have Ach > 100. After the high Ach requirement, the atmospheric neutrino 
simulation peaked at 10 TeV, while the signal simulation peaked around 100 TeV (Figure 
The energy range defined by the central 90% of the signal with A^^h > 100 is the energy 
range for the sensitivity or limit. For this search, the central 90% signal region extends from 
16 TeV to 2.5 PeV. 

The efficiency of the detector for neutrinos is quantified by the effective area. In the 
energy range relevant to this analysis, it increases with energy and is further enhanced by 
including uncontained events. The effective area is described by the following equation where 
A^ represents the number of observed events and T is the detector livetime: 



13 



Data 2000 - 2003 

Uncertainty in atms. v 

Barr et al. atms. v 

Honda et al. atms. v 

CORSIKAatms.n 

E^d) = 1 0"^ GeV cm-^ s'^ sr-^ 

sig 



« 

£|io« 

10^ 
10^ 
10 



wlO^ 



^ 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 


L level 2 


• 

• :"!L 











-0.8 -0.6 -0.4 -0.2 -0 
Cosine of Reconstructed Zenith Angle 





-1 -0.8 -0.6 -0.4 -0.2 -0 

Cosine of Reconstructed Zenith Angle 



«10^ 

c 

* 

£iio« 

10* 
10" 
10^ 
10=^ 
10 



level 3 



-0.8 -0.6 -0.4 -0.2 -0 
Cosine of Reconstructed Zenith Angle 



-1 -0.8 -0.6 -0.4 -0.2 -0 

Cosine of Reconstructed Zenith Angle 




-1 -0.8 -0.6 -0.4 -0.2 -0 

Cosine of Reconstructed Zenith Angle 



FIG. 3: The cosine of the zenith angle is plotted for all events surviving the event quality criteria 
at a given level. Events at cos(zenith) = — 1 are traveling straight up through the detector from 
the Northern Hemisphere. The initial zenith angle requirement removed events from 0° to 80" 
(level 1 - top right). Events reconstructed just above the horizon appear at the right side of each 
plot. Each level represents an increasingly tighter set of quality requirements. As the quality level 
increased, misreconstructed downgoing muons were eliminated. To ensure a clean upgoing sample, 
events coming from the horizon were discarded by requiring reconstruction angles greater than 
100°. The final analysis was performed at level 5 (bottom right) with horizontal events removed. 
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FIG. 4: Nch, or number of OMs hit. Prediction for both conventional and prompt atmospheric 
neutrinos are shown and their uncertainties are represented by the gray band. The central prompt 
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32( 1 models. 



neutrino flux used here is the average of the Martin GBW [30] and Naumov RQPM 
All atmospheric neutrinos are normalized to the number of data events in the range 50 <A'ch <100. 
The lower signal flux curve corresponds to the limit obtained in this paper. 




FIG. 5: The true energy of the simulation is shown for atmospheric neutrino and signal events. 
The thin dashed (atms. v) and solid (signal ly) curves represent the number of events before the 
-^ch > 100 requirement. The thick dashed and solid lines represent only the events in the high 
energy sample. 
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FIG. 6: Effective area for as a function of the true simulated energy at the Earth's surface in 
intervals of cosine of the true zenith angle. The effective area is the equivalent area over which 
the detector would be 100% efficient for detecting neutrinos. The absorption of neutrinos and 
reduction of their energy via neutral current interactions in the Earth are taken into account. The 
angle-averaged effective area is represented by the solid black line. 



f = / A''^Q{E,,n)^,dndE . (1) 

The effective BjrGcl clS db function of energy is shown for different zenith angle regions in 
Figure E] (and is tabulated in Appendix B). At energies greater than 10^ GeV, the Earth 
begins to be opaque to neutrinos depending on direction and the highest energy events are 
most likely to come from the region around the horizon 28||. In Figure El the effective 



area decreases at high energy because tracks with zenith angles between 80° and 100° were 
discarded. Most of the events that were removed were high energy events from the horizon. 



V. SYSTEMATIC UNCERTAINTIES 



A discovery is made if an excess of events over the predicted background is observed in the 
data. However, due to uncertainties in the simulation, the number of signal and background 
events predicted may not accurately reflect the true signal and background. Theoretical 
uncertainties exist in the atmospheric neutrino flux models for several reasons. The cosmic 
ray spectrum is very uncertain at high energy and hadronic interactions in this energy range 
are not well understood. There are also detector-related uncertainties. Photons scatter more 
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in dirty or bubble-laden ice. Hence, our incomplete understanding of the dust layers in the 
ice and the bubbles in the hole ice (formed from water that refroze after deployment of the 
OMs) add uncertainty to our models 29|]. There are also uncertainties in the simulation 
associated with the modeling of light propagation in the ice and with the optical module 
sensitivity. These contributions are considered individually to see how they affect the number 
of simulated events in the final sample. The number of experimental data events remaining 
after the final energy requirement (A^^ch > 100) is then compared to the range of predicted 
background and signal events when uncertainties are considered. 



A. Theoretical Uncertainty in the Background 



For this analysis, two models based on the work of Barr et al. [21|, |22|, |23| and Honda 
et al. were considered equally likely options for the background atmospheric neutrino 
simulation. These two models are recent calculations that cover the highest and lowest 
portion of the atmospheric neutrino flux band created by uncertainties in the primary cosmic 
ray flux and the high energy hadronic interaction models. Since these models do not extend 
to the high energies needed for this analysis, the models were extrapolated to higher energies 
based on the procedure described in Appendix C. Differences between the Barr et al. and 
Honda et al. models are also summarized there. 

Conventional atmospheric neutrinos from the decay of pions and kaons are not the only 
source of atmospheric background. Above 50 TeV - 1 PeV, the source of atmospheric 



neutrinos is expected to change 



30 



31 



32 



33, 



34| . Semileptonic decays of short-lived 



charmed particles become the main contributor to the atmospheric neutrino flux. Since these 
charmed particles decay quickly before they can lose much energy, the resulting neutrinos are 
called prompt neutrinos. At these energies, charm quarks are produced primarily via gluon- 
gluon fusion. Uncertainties in the gluon distribution at low Bjorken x lead to uncertainties 
in this prompt lepton flux. 

Uncertainties were included for both conventional atmospheric neutrino models. The 
uncertainty in the cosmic ray spectrum was estimated as a function of energy based on 
the spread of values measured by many cosmic ray experiments 35|]. These uncertainties 



were added in quadrature with the estimated uncertainty due to choosing different hadronic 



interaction models l2ll . l36l |. Uncertainties were also estimated based on the spread of 
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predictions surrounding the unknown prompt neutrino flux. Unless mentioned otherwise, 
when prompt neutrinos were included in this work, the average of the Martin GBW (Golec- 
Biernat and Wiisthoff) 30| and Naumov RQPM (Recombination Quark Parton Model) 
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32| models is shown. This is henceforth called the central prompt neutrino model. 



All of the uncertainty factors for the total (conventional + prompt) atmospheric neutrino 
flux were combined and are shown as a function of energy in Appendix C. Since the true 
energy of every simulated event is known, each event was given a weight based on the 
maximum uncertainty estimated for that neutrino energy. As a result, three predictions 
for the number of atmospheric neutrinos in the final high energy sample were made (the 
model, the model plus maximum energy-dependent uncertainty, the model minus maximum 
energy-dependent uncertainty). Since both the Barr et al. and Honda et al. fluxes were 
considered equally likely, the central prompt neutrino flux was added to both predictions. 
Then uncertainties were added and subtracted to both of these total atmospheric neutrino 
fluxes, creating six different background possibilities. 



1. Normalizing the Atmospheric Neutrino Simulation to the Data 

After all but the A^^h event selection requirements were fulfilled, the A^h distribution for 
the observed low energy events was inconsistent with that for the total atmospheric neutrino 
simulation in normalization. Each of the six atmospheric neutrino background predictions 
was renormalized to match the number of data events observed in the low A^h region, where 
the signal was insignificant compared to the background. By rescaling the simulation to 
the number of observed data events, the uncertainty in the atmospheric neutrino flux was 
reduced to the uncertainty in the spectral shape. 

Since some of the atmospheric neutrino models predicted more events than the data while 
others predicted less, renormalization of the models to the data brought the simulated models 
into closer agreement. The renormalization is explained in greater detail in Appendix C. 

Since the purpose of this normalization was to correct for theoretical uncertainties in the 
atmospheric neutrino background prediction, it was not applied to the simulated neutrino 
signal. 
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FIG. 7: In the inverted analysis, the highest quahty downgoing events were studied. The A^ch 
distribution is shown for all events which survive the inverted quality requirements. 

B. Simulation Uncertainties 

To assure that the detector response to high energy events (A^ch > 100) is understood, it 
is important to study energetic events while simultaneously keeping the high energy upgoing 
events blind to the analyzer. To this end, an inverted analysis was performed in which high 
quality downgoing tracks were selected from the initial data set. The advantage of studying 
high quality downgoing tracks is that large data sets are available to study both the high and 
low energy events. When the data and simulation observable distributions are not perfectly 
matched, imposing event quality requirements may result in removing different fractions of 
the simulation in comparison with the data. The inverted analysis was used to study this 
systematic effect. 



1. Inverted Analysis Using Atmospheric Muons 

For the inverted analysis, all event quality requirements described previously (and sum- 
marized in Appendix A) were applied, but events were selected based on a high probability 
of being downgoing rather than upgoing tracks. 

When compared to the downgoing experimental data, small shifts were observed in the 
peaks of the simulated distributions for the number of direct hits (A^dir), the smooth distribu- 
tion of hits along the track (S"), the event-by-event track resolution {MR) and likelihood of 
being downgoing muon tracks rather than upgoing (Inverted Likelihood Ratio, ILR). These 
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discrepancies are most likely due to inaccurate modeling of optical ice properties in the 
simulation, since it is technically challenging to implement a detailed description of photon 
propagation through layered ice. 

If multiple parameters are correlated, it is possible that mismatches in one parameter 
may affect the agreement between data and simulation in another. In order to study these 
effects, the differences in the data and simulation were analyzed at the level where no quality 
criteria had been applied. The simulated distributions needed to be shifted to larger values 
by approximately 10% for A^dir, 8% for S, 5% for MR and 1% for ILR. When simultaneous 
corrections to the simulation for all of these effects were applied, the downgoing data and 
simulation were in better agreement for all parameter distributions. Later in the analysis, 
these shifts were applied to the upgoing analysis. The number of background and signal 
events appearing in the final upgoing sample was recalculated based on these simulation 
shifts. 

2. Uncertainties in Detector Response 

The downgoing sample from the inverted analysis was also used to study how well the 
detector response was simulated in the high energy {N^h > 100) regime. Using downgoing 
data and atmospheric muon simulation, a ratio of the number of events was taken as a 
function of A'^ch from the histograms shown in Figure [71 If the simulation perfectly described 
the data, the shapes of the N^h distributions would match and this ratio would be flat. The 
downgoing ratio was mostly fiat, but slightly increased at large N^^ where low statistics 
introduced large uncertainties. The statistical uncertainty aside, a scenario was considered 
in which the downgoing data to simulation ratio truly increased as A^ch increased. Under 
this scenario, the simulation is renormalized by a larger factor at high A^^ch to replicate the 
data. This A"ch-dependent renormalization was then applied to the upgoing simulation used 
for the main part of the analysis. This non-linear normalization factor had a negligible 
effect in the number of atmospheric neutrinos predicted in the final sample of events with 
A"ch ^ 100. However, the high energy signal simulation event rate increased by 25% when 
this non-linear A"ch effect was included. This uncertainty was incorporated in the final limit 
calculation that will be described in the next section. 

Detection efficiency also depends on the OM sensitivity. This parameter of the simulation 
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was modified and new simulated events were generated. After comparing tlie data and 
simulation with different OM sensitivities, a 10% uncertainty in the total number of events 
due to inaccurate modelling of the OM detection sensitivity was incorporated into the final 
upper limit calculation. 

The systematic errors due to the neutrino interaction cross-section, rock density (below 
the detector), and muon energy loss do not contribute significantly to this analysis |16i] . 



VI. RELATIONSHIP BETWEEN UPGOING AND DOWNGOING EVENTS 



In addition to using the inverted analysis to study high energy events and the bias in- 
troduced by inaccurate simulation, the downgoing events can be used as a calibration beam 
for the upgoing atmospheric neutrino flux. To do this, the same simulation package (COR- 
SIKA v6.030, QGSJETOl, ^primary oc E~^-'^) was used to describe the downgoing atmospheric 
muons and the upgoing atmospheric neutrinos [37 1. 

As shown in Table [H the ratio of experimental data to CORSIKA downgoing muon sim- 
ulation was relatively constant as the event selection became more discriminating. The 
simulation does not match the data normalization and this may be a consequence of the 
theoretical imperfections in the CORSIKA simulation (mainly due to the hadronic interac- 
tion model (QGSJETOl) and uncertainty in the primary spectrum ($ oc E"^'"^)). Another 
contributing factor to the normalization difference may be that light propagation in the 
layered ice is modeled inaccurately. When the upgoing CORSIKA atmospheric neutrinos 
are rescaled based on the downgoing muons, then the upgoing experimental data and COR- 
SIKA atmospheric neutrino simulation are in good agreement for the number of low energy 
events in the final sample. This can only be seen when the tightest criteria are applied be- 
cause misreconstructed muons and coincident muons contaminate the data sample when the 
quality requirements are loose. For instance, at level 5 in the inverted analysis, the ratio of 
downgoing data to simulation was 1.22. For the upgoing analysis at level 5, 146 events were 
observed and 124.9 CORSIKA atmospheric neutrinos were predicted. When adjusted based 
on the inverted analysis, 152.4 (= 124.9 x 1.22) CORSIKA atmospheric neutrinos were pre- 
dicted, which is in good agreement with the observed value. This shows that it is possible to 
adjust the normalization of the upgoing events based on the downgoing observations (when 
the up and downgoing simulation use the same input assumptions about the spectrum and 
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LI 


L2 


L3 


L4 


L5* 


Downgoing 












data (xlO^) 


7.88 


6.70 


6.05 


5.89 


2.59 


CORSIKA 












alms. //(xUf) 


(3.63 


5.75 


5.12 


5.01 


2.12 


ratio 


1.19 


1.17 


1.18 


1.18 


1.22 


Upgoing 












signal 


325 


241 


191 


185 


103 


coinc /J, 


2570 


268 


45.8 


29.4 





misreconstructed 












CORSIKA 












atms. fi 


37800 


2570 


148 


34.2 





Barr et al. 












atms. u 


681 


526 


393 


380 


194 


Honda et al. 












atms. v 


513 


400 


300 


290 


149 


Martin GBW prompt v 


1.9 


1.9 


1.6 


1.5 


0.7 


Naumov RQPM prompt v 


18.9 


18.9 


16.0 


15.5 


7.5 


CORSIKA 












atms. u 


440 


335 


251 


243 


125 


Adjusted 












CORSIKA 












atms. v 


524 


392 


296 


286 


152 


data 


276894 


24422 


1269 


531 


146 



*L5 = level of final analysis 



TABLE I: The number of low energy events (50 <A'ch <100) at a given level (see Appendix A) 
for the different types of simulation and experimental data. The top portion of the table presents 
results from the inverted analysis. The main upgoing analysis is summarized in the lower portion 
of the table. Note that the upgoing data and adjusted CORSIKA atmospheric neutrino flux are in 
good agreement when the CORSIKA neutrino events are adjusted by the scale factor determined 
in the downgoing analysis. This agreement can be seen at the tightest quality levels because all 
misreconstructed backgrounds have been removed. 

interaction model). 



VII. RESULTS 

We calculated a confidence interval based on the number of events in the final AT^h > 100 
sample of the predicted background and signal and the observed data. Statistical and 
systematic uncertainties were incorporated into the confidence interval such that the true, 

22 



Systematic Uncertainty 


Theoretical uncertainty in atms. v flux 


See Figure [12] 


Number of Direct Hits 


10% 


Smoothness 


8% 


Median Resolution 


5% 


Inverted Likelihood Ratio 


1% 


Total background uncertainty 


+19% / -18% 


Non-linear detector response at high A^ch 


25% 


OM sensitivity 


10% 


Total signal efficiency uncertainty 


+/- 27% 



TABLE II: The systematic error was estimated with several techniques. The theoretical uncertainty 
in the atmospheric neutrino flux was estimated as a function of energy (Section IV Ap . Using the 
inverted analysis, shifts were observed between the data and simulation in four parameters (Section 
IVB ip . When each of the above mentioned uncertainty factors was applied to the atmospheric 
neutrino simulation, the resulting spread in the number of events predicted in the Nf.\i > 100 
sample indicated that the total background uncertainty was +19% / -18% of the average predicted 
background, 7.0 events. The non-linear response of the detector in Net was estimated as 25% 
(Section IVB2p . When added in quadrature with the 10% uncertainty in OM sensitivity (Section 
IVB 2p . the total signal efficiency uncertainty was +/-27%. 



but unknown, value of the diffuse flux of astrophysical neutrinos is contained within the 
interval in 90% of repeated experiments. A hybrid frequentist-Bayesian method based on the 
work of Cousins and Highland |S] was used to construct a confidence belt with systematic 
uncertainties. The likelihood ratio ordering was based on the unified confidence intervals 
explained by Feldman and Cousins 27(1 . The uncertainty in the detection efficiency of the 
signal was set at 27% (10% for optical module sensitivity added in quadrature with 25% 
for non-linearity in the A'ch spectrum when data and simulation are compared). Systematic 
uncertainties on the number of background events in the final sample were also included 
in the confidence belt construction. Inclusion of the signal and background uncertainties 
followed the methods described by Conrad et al. [39] and Hill 40] . 



In constructing the fiat Bayesian prior for the background, twelve atmospheric neutrino 
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models were considered equally likely. The twelve predictions were derived as follows. Ini- 
tially, two background predictions were considered, Barr et al. and Honda et ai, each with 
the central prompt neutrino flux added. To include systematic uncertainties in the models, 
maximum uncertainties were added and subtracted from each model (Section IV Ap . Hence, 
the six predictions were named Barr et al. maximum, nominal and minimum and Honda et 
al. maximum, nominal and minimum. The number of events predicted for the background 
in the flnal sample is listed in Appendix C. To account for systematic uncertainties in the 
detector response, the simulation was shifted in four different parameters as described in 
Section IVB 1[ This simulation shift was performed on each of the 6 models described above, 
hence creating a total of 12 different atmospheric neutrino predictions that were used in the 
confldence belt construction. The number of events predicted by the 6 models with shifted 
simulation was within 10% of each number reported in Appendix C. 

A. Results for $ cx 

The signal hypothesis consisted of a flux = 1.0 x 10^^ GeV cm~^ s~^ sr^^. At this 
signal strength, 66.7 signal events were expected in the flnal A^ch > 100 data. (This value 
assumes half of the correction from the simulation shifts since 68.4 events were predicted in 
the flnal selection, but the number of events decreased to 65.0 when the simulation shifts were 
applied.) The sensitivity was obtained from the slice of the confldence belt corresponding to 
zero signal strength. The median observation assuming no signal was seven events, giving a 
median event upper limit of 6.36 and hence a sensitivity of 9.5 x 10^® GeV cm~^ s^^ sr^^. 

When the data with Nch > 100 were revealed, six data events were observed. This was 
consistent with the average expected atmospheric neutrino background of 7.0 events (after 
averaging all models that have been rescaled to the low energy data). Information about 
the observable quantities for the flnal six events can be seen in Table IIIII The flnal Net 
distribution is shown in Figure |H The total number of events predicted for the signal and 
background can be compared to the observed data in Table |Vl] {Nch <100) and in Table 
IVIII {Nch > 100). With uncertainties included, the upper limit on a diffuse $ oc flux of 
muon neutrinos at Earth (90% confldence level) with the AMANDA-II detector for 2000 - 
2003 is 7.4 X 10^^ GeV cm~^ s^^ sr^^ for 16 TeV to 2.5 PeV. The results are compared to 
other neutrino limits in Figure [HI 
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FIG. 8: The upper limits on the flux from sources with an E energy spectrum are shown for 
single and all-flavor analyses. All-flavor upper limits have been divided by three, assuming that 
the neutrino flavor ratio is 1:1:1 at Earth. The Frejus [12], MACRO [13[, and AMANDA-BIO 



lai 

3 



upper limits on th e u, ,, flux are shown, as well as the unfolded atmospheric spectrum from 2000 



AMANDA-II data 



43]. The AMANDA-II all-flavor limit from 2000 11 1, the AMANDA-BIO UHE 



3 



and the RICE six year limit 



44| have all been adjusted for 



limit [43] , the Baikal flve year limit 
the single flavor plot. The ^ oc E^^ limit from this analysis is a factor of four above the Waxman- 
Bahcall upper bound. Although not shown, this analysis excludes the oc E~^ predictions made 
by Nellen, Mannheim and Biermann [4] and Becker, Biermann and Rhode 5| and constrains the 
MPR upper bound for optically thick pion photoproduction sources Gj. The IceCube sensitivity 



for a full detector was estimated with AMANDA software 



45]. 



B. Results for Other Energy Spectra 



Other signal models were also tested with this data set. Due to their different energy 
spectra, the A^ch requirement was reoptimized by minimizing the MRF with each signal 
model. For signal models with softer spectra than $ oc E~^, a lower Net requirement was 
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Event 


1 


2 


3 


4 


5 


6 


required 
value 


Year 


2001 


2001 


2001 


2001 


2002 


2003 




Day of Year 


118 


186 


210 


274 


226 


182 






102 


106 


157 


116 


100 


111 


>100 


Track Length [m] 


206.7 


221.8 


197.7 


178.2 


180.4 


207.6 


>170 


Number of 
Direct Hits 


27 


32 


30 


22 


29 


29 


>13 


Zenith Angle [°] 


107.3 


121.6 


106.1 


101.8 


123.8 


113.3 


>100 


Median Resolution [°] 


2.4 


1.4 


1.8 


3.0 


1.6 


2.8 


<4.0 



TABLE III: Observable and reconstructed qualities are shown for the final six events. In addition, 
events fulfilled requirements based on the reconstructed values of their smoothness (S) and their 
upgoing vs. downgoing likelihood ratios. 



Experiment 


Upper Limit 

[GeV cm^^ s^^ sr^^] 


Energy Range 

logio [K (GeV)] 


Muon neutrinos only 


Frejus [12] 


5.0 X 10"^ 


~3.4 


MACRO [13] 


4.1 ±0.4 X 10"^ 


4.0 - 6.0 


AMANDA-BIO [10] 


8.4 X 10"' 


3.8 - 6.0 


AMANDA-II (this analysis) 


7.4 X 10"*^ 


4.2 - 6.4 


All neutrino flavors 


Baikal [14] 


8.1 X lO"'^ 


4.3 - 7.7 


AMANDA-BIO [43] 


0.99 X 10"^ 


6.0 - 9.5 


AMANDA-II [llj 


8.6 X 10"' 


4.7-6.7 



TABLE IV: Upper limits for the diffuse flux of extraterrestrial neutrinos as reported by a number 
of experiments. The first four analyses only constrain the flux of v^, + , while the last three 
constrain the total neutrino flux, (z^e + i^e + ^^/^ + ^/x + + i^r)- 



optimal, > 71. Four prompt neutrino models 30|, l3l|, l32, l3J] and one astrophysical neu- 
trino model \^ were tested under these conditions. One astrophysical model was optimized 
at A^ch > 86 [6]. Two astrophysical neutrino models with harder spectra than $ oc E"^ were 

n p n 

tested with a higher energy requirement, A^ch > 139 [6|, 17|, 18|. 

Results of these searches are summarized in Table |Vl The normalization of the over- 
all number of low energy atmospheric neutrinos to data was performed over the re- 
gion 50 < A'ch < 100 for the harder spectra (A^^h > 139), and over 50 < A^ch < 71 and 
50 < A^ch < 86 for the softer spectra. 

When the data from the A^ch > 139 region were examined, there was good agreement 
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with the expected atmospheric neutrino background (1 event observed on a backround of 
1.55). For A^^ch > 86, 14 events were observed while an average of 12.9 background events 
were predicted. However, 37 events were observed while only 27.4 events were expected 
for A'ch > 71, leading to a two-sided confidence interval. Since the chance probability of 
observing 37 or more events on this background is 4%, we do not exclude the background- 
only null hypothesis. The 90% confidence interval for fi is shown for each model in Table 
IVl and upper limits are calculated based on the upper bound of each confidence interval. If 
the MRF is greater than 1, then the model is not ruled out based on observations from this 
four- year data set. Since more events were observed in the data than were predicted by the 
background simulation for A^ch > 71, the upper limit on those five models is roughly a factor 
of three worse than the sensitivity. 



1. Astrophysical Neutrinos 

The first astrophysical neutrino model tested with the A'^ch > 139 requirement was initially 

proposed by Stecker, Done, Salamon and Sommers [tI. The flux tested in this analysis 

includes the revision in the erratum of their original paper Ml and the factor of 20 reduction 
fl 

by Stecker in 2005 [8||. The model predicts a flux ($sdss) of high energy neutrinos from the 
cores of AGNs, especially Seyfert galaxies. Based on the present data, the upper limit on 
this flux is 1.6-$sDss- The best previous limit on this model was established by the Baikal 



experiment, with an upper limit of 2.5 -^sdss 

Mannheim, Protheroe and Rachen (MPR) 6| computed an upper bound for neutrinos 
from generic optically thin pion photoproduction sources (r„^ < 1), as well as an upper 
bound for neutrinos from AGN jets. (In addition, they calculated an upper bound for generic 
optically thick (t„^ ^ 1) pion photoproduction sources assuming a $ oc E^^ spectrum, but 
this is constrained by the results discussed in the previous section.) The upper bounds do 
not necessarily represent physical neutrino energy spectra, but were constructed by taking 
the envelope of the ensemble of predictions for smaller energy ranges. Each flux prediction 
within the ensemble was normalized to the observed cosmic ray proton spectrum. 

Nonetheless, the shapes of these two upper bounds were tested as if they were models. 
However, one should be careful not to misinterpret the results. A limit on a model implies 
a change in the normalization of the entire model. A limit on an upper bound only implies 
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a change in normalization of the bound in the energy region where the detector energy 
response to that spectral shape peaks. 

The MPR AGN jet upper bound was tested with the A'^^h > 139 requirement. The upper 
limit on this spectrum is 2.0-$mpragn- In comparison, the Baikal upper limit on this 
spectrum is 4.0-<I>mpragn- 

The MPR upper bound for optically thin sources was tested with a Nch > 86 requirement. 
The limit on this spectrum and normalization is 0.22-$MPRr<i- 

The remaining neutrino searches were conducted with the lower A/'ch^requirement, 
Nch > 71. A signal hypothesis involving neutrinos from starburst galaxies [9'] was tested. 
Loeb and Waxman assumed that protons in starburst galaxies with energy less than 3 PeV 
convert almost all of their energy into pions. Their work predicts a range that should en- 
compass the true neutrino spectrum, but the model tested here uses the most probable 
spectrum from the paper, $ oc E~^'^^. This analysis assumed the flux was valid for energies 
ranging from 10'^ to 10^ GeV. The upper limit on this spectral shape and normalization is 

21.1-$starburst- 

These astrophysical neutrino models and their observed upper limits based on this data 
set are shown in Figure [91 Neutrino oscillations are taken into account for all models where 
this factor was not already applied. 

2. Prompt Neutrinos 

Since prompt neutrinos have a harder (less steep) spectrum than the conventional atmo- 
spheric neutrinos, it is possible to search for a prompt neutrino flux by separating the two 
event classes in energy. The final A^h requirement was reoptimized yielding A^h > 71 and 
the normalization factor was determined based on the interval (50 < Ach < 71). 

In the astrophysical neutrino searches described thus far, the range of atmospheric neu- 
trinos predicted in the final sample included an uncertainty due to the unknown prompt 
neutrino flux. For the search for prompt neutrinos, this uncertainty in the total atmospheric 
neutrino flux was changed so that only conventional atmospheric neutrino uncertainties were 
included. Since the atmospheric neutrino simulation was still normalized to the low energy 
data, the overall effect in the atmospheric background prediction for the final sample was 
small. 
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FIG. 9: Astrophysical neutrino models and upper limits established with this analysis. The Barr 
et al. and Honda et al. atmospheric neutrino models are shown as thin lines with maximum 
uncertainties assumed by this analysis represented by the band. Other models that were tested 



includec 



sources 



the SDSS AGN core model 



re MPR upper bounds for AGN jets and optically thin 



6| , and a starburst galaxy model 



a. 



Martin et al. predict prompt lepton fluxes based on the GBW model for deep inelastic 
scattering. This model includes gluon saturation effects 30|] which lower the predicted charm 
production cross sections. The predicted flux is lower than the sensitivity of this data set. 
The upper limit on this model is 60.3 • $MartinGBW- 



The Naumov RQPM 
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32l | model of prompt atmospheric neutrinos incorporates data 



from primary cosmic ray and hadronic interaction experiments. This non-perturbative model 
includes intrinsic charm [33|. The upper limit on this model is 5.2-$NaumovRQPM- 



Prompt neutrinos based on the models of Zas, Halzen and Vazquez were also simulated 



3J]. A parameterization was established to describe the energy dependence of the charm 



cross section. For the Charm C model, the charm cross section was fltted to experimental 
data. In the Charm D model, the cross section was parameterized by Volkova [4l|. Due 
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FIG. 10: Prompt neutrino models and upper limits based on this analysis. The Barr et al. and 
Honda et al. atmospheric neutrino predictions are shown for reference. Two charm models 
were tested, along with the Naumov RQPM 
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32] and Martin GBW 



30l | models. 



to the upward fluctuation in the number of events in the A^ch > 71 region, the upper hmit 
for Charm C is l.S-'I'charmC- The upper hmit on the Charm D model is O.OS-'I'charmD- The 
MRF is less than 1.0, hence the Charm D model is disfavored at the 90% confldence level. 

The prompt neutrino models are shown in Figure [ini along with the upper limits based 
on these data. 



VIII. CONCLUSIONS 



The experimental data were consistent with the predicted range of atmospheric neutrino 
background. Six high energy events were observed in the flnal data set, while the average 
predicted background was 7.0 events. There is no indication of an astrophysical signal. 
At a 90% confidence level, the diffuse flux of extraterrestrial muon neutrinos with an 
spectrum is not larger than 7.4 x 10^^ GeV cm~^ s~^ sr~^ for 16 TeV - 2.5 PeV. 

This analysis also provides upper limits on four astrophysical neutrino models and four 
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upper limit 
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0.074 x$E-2 


1.6 X$SDSS 


2.0 X$MPRAGN 


{logioEmin, logioEjnax) 
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(5.0,6.9) 





MPR Tn-y < 1 \6] 
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^obs 


14 


37 
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(0,9.49) 


(0,22.13) 
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h/hs X $ 


0.22 X$MPRr<l 


21.1 X <i>Starburst 


{logioEmin, logioEmax) 


(4.0,5.8) 


(3.8,6.1) 



Prompt u 
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Naumov 








GBW [30] 


RQPM [31, 32] 


CharmC [34] 


CharmD [34] 


iVeh 


71 


71 


71 


71 


rih 


27.4 


27.4 


27.4 


27.4 


Us 


0.41 


4.74 


16.05 


26.15 


/^median (^b ) 


8.75 


8.75 


8.75 


8.75 


sensitivity 










Mmedian(nb)/ns X $ 


21.3 X$MGBW 


1.8 X$NRQPM 


0.55 X<i)charmC 


0.33 X$CharmD 


^obs 


37 


37 


37 


37 


/^90%C.I. 


(1.29,24.72) 


(1.29,24.72) 


(1.29,24.72) 


(1.29,24.72) 


upper limit 










^i/Us X $ 


60.3 X$MGBW 


5.2 X$NRQPM 


1.5 X$ChannC 


0.95 X$CharmD 


{logioEmin, logioEmax) 


(3.5,5.5) 


(3.6,5.6) 


(3.8,5.7) 


(3.6,5.6) 



TABLE V: Several flux shapes were tested with this data set. N^h is the minimum number of OMs 
that had to be hit for an event to appear in the final data set. The predicted number of events 
for background, rib, and signal, rig, were determined by the simulation. The median event upper 
limit is /Umedian(f^b)- The Sensitivity is the model flux multiplied by the median event upper limit 
and divided by the number of signal predicted. The number of events observed in the four year 
data sample is nobs- The upper limit is calculated from the maximum value of the 90% confidence 
interval for the event upper limit, fi. The upper limit is the test flux multiplied by fi/ug- All values 
quoted here incorporate systematic uncertainties. 
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prompt neutrino models. For the hardest signal spectra, the results are consistent with 
background. The softer spectra were tested with lower A'ch requirements and despite the ob- 
servation leading to a two-sided 90% confidence interval, the level of excess is not significant 
enough to claim a detection. 

Before requiring events to fulfill N^h > 100, the observed events were compared to the 
atmospheric neutrino simulation with systematic uncertainties included. The observed low 
energy data were used to normalize the atmospheric neutrino simulation, hence narrowing 
the range of atmospheric neutrinos predicted by the different models for the final high energy 
sample. Systematic effects of the event selection procedure were studied in the inverted 
analysis using atmospheric muons. A consistency was established between the observed 
downgoing atmospheric muon flux and the upgoing atmospheric neutrino flux using the 
inverted analysis. 

This result is the best upper limit on the diffuse flux of muon neutrinos to date. The 
upper limit is an order of magnitude lower than the previous AMANDA result by performing 
a multi-year analysis 10| and by using a larger detector, AMANDA-II instead of AMANDA- 
BIO. For a $ oc E^^ spectral shape, this analysis provides an upper limit that is a factor 
of three better than the Baikal muon neutrino upper limit (muon neutrino upper limit = 
all-flavor limit /3 assuming a 1:1:1 flavor ratio). 

This analysis set upper limits on four prompt atmospheric neutrino predictions, while one 
of these models is disfavored at a 90% confldence level. Other spectral shapes were tested 
for astrophysical neutrinos. No models were excluded, however constraints were placed on 
the existing predictions. The shapes of the MPR upper bounds were tested in the energy 
region where the detector response peaks. For the benchmark $ oc E^^ spectral shape, the 
current limit is a factor of 4 above the Waxman-Bahcall upper bound. 

AMANDA-II has now been integrated into IceCube. The sensitivity of the IceCube 
detector will continue to improve as the detector grows to its flnal volume, 1 km^. Based on 
estimations with AMANDA software, the full IceCube detector will have a sensitivity that 
is a factor of 10 better than this analysis after one year of operation 45 1. 
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Upgoing <iVch <100 




LO 


LI 


L2 


L3 


L4 


L5 


signal 
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coinc. /X 




60800 
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misreconstructed 
CORSIKA 

atms. /i 
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Barr et al. 
atms. v 
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4470 


3890 
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Honda et al. 

atms. v 
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3130 


420 


Martin GBW 
prompt atms. v 
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8.2 


6.4 
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1.2 


Naumov RQPM 
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71.6 


56.7 


50.7 


11.5 


Data 
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4088 
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Downgoing <Nch <100 


CORSIKA 
atms. jj, (xlO'') 


386 




288 


212 


195 


24 


data (xlO'^) 


432 




323 


255 


229 


30 



L5 = level of final analysis 



TABLE VL The number of low energy events (0 <Nch <100) at a given quality level for the 
different types of simulation and experimental data. 
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Upgoing A^ch > 100 
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Naumov RQPM 
prompt atms. v 
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4.8 


4.2 
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2.2 


Data 




11456 


1347 


96 


45 
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Downgoing iVch > 100 


CORSIKA 
atms. II (xlO^) 


7.31 




6.53 


6.05 


6.01 


5.09 


data (xlO'^) 


9.75 




8.59 


8.07 


8.03 


6.60 



L5 = level of final analysis 



TABLE VIL The number of high energy events {N^h >100) at a given quality level for the different 
types of simulation and experimental data. 
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Appendix A: Event Selection Techniques 

Event selection techniques were applied to find the best reconstructed upgoing tracks. 
The event requirements were tightened through a series of values, becoming more restric- 
tive at each of the five different levels. As seen in Figure [TT], requiring a minimum value 
of the track length, for instance, can be a powerful method of rejecting misreconstructed 
downgoing backgrounds. The event selection requirements for L^i^, A^dir? smoothness, me- 
dian resolution and likelihood ratio were established to remove many orders of magnitude 
more misreconstructed background than upgoing atmospheric neutrinos or signal neutrinos. 
Events which did not meet an optimized minimum or maximum value of each parameter 
were removed. 

The strength of these quality requirements was adjusted at each level. The requirement 
is defined for each parameter in Table IVIIII The plots in Figure |3] show the zenith angle 
distribution of all events fulfilling the zenith angle >80° and event observable requirements 
at the chosen level. After the zenith angle criteria was fulfilled at Level 1, the data mostly 
contains misreconstructed atmospheric muons (top right. Figure [3]) • As the quality parame- 
ters become more restrictive, the data begins to follow the atmospheric neutrino simulation 
in the upgoing direction and the atmospheric muon simulation in the downgoing direction. 
At Level 5, the event quality requirements were strong enough to have removed all of the 
misreconstructed downgoing atmospheric muon events that were simulated. However, to be 
sure that the final data set only included atmospheric and astrophysical neutrinos and no 
misreconstructed downgoing events, an additional zenith angle requirement was imposed. 
All events were kept if they were reconstructed between 100° and 180°. The analysis con- 
tinued with the data sample shown at Level 5. 
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FIG. 11: The reconstructed track length within the detector is shown. In order to identify muon 
neutrino tracks, events were required to have long tracks of at least 170 meters. This removed a 
large fraction of the atmospheric muon simulation, but had a smaller effect on the atmospheric 
neutrino and signal simulations. 
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<0.30 
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Resolution [°] 










<4.0 


<4.0 
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AL > -38.2 cos{zenith) + 27.506 


Number of 
Remaining Events 


5.2 X 10'' 


7.8 X lO'' 


1.2 X 10"^ 


3.5 X 10^ 


1.8 X lO'"^ 


465 



* = level of the final analysis 



TABLE VIII: The table summarizes the event quality requirements as a function of quality level. 
Events only remained in the sample if they fulfilled all of the parameter requirements for a given 
level. The removal of all horizontal events (zenith < 100) contributed to the large decrease in 
events from L4 to L5. 

Appendix B: Neutrino Effective Area 
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TABLE IX: Effective area as a function of the energy and zenith angle of the simulation. 
Appendix C: Atmospheric Neutrino Flux 

For this analysis, the atmospheric neutrino flux models by Barr et al. and Honda et 
al. were both considered equally likely options for the background atmospheric neutrino 
simulation. These two models are among many that use slightly different parameterizations 
of the all-nucleon cosmic ray flux to derive the atmospheric neutrino flux. 
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49 


7.8 


35.5 
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TABLE X: The angle-averaged neutrino effective area as a function of energy. 



For this analysis, the Barr et al. flux below 10 GeV was taken from 211]. From 10 GeV 



to 10 TeV, the flux tables from 



23|, were used. 



221 ]. based on the primary spectrum of 
Above 10 TeV, the weight was derived by performing a 2-dimensional flt with a flfth degree 
polynomial to the logioE vs. cos(zenith) tables of the atmospheric neutrino flux values from 
lower energies just mentioned. The TARGET version 2.1 [46| hadronic interaction model 
was used |2l| . 

In an attempt to better flt the AMS [l^] and BESS Q, Q data, Honda et al. changed 
the power law flt to the proton cosmic ray spectrum from -2.74 to -2.71 above 100 GeV 20]. 
Other parameters in the cosmic ray flt remained similar to the Barr et al. flux mentioned 
above 36|], although the DPMJET-III 50|] interaction model was used. The atmospheric 



neutrino weights from [20] were used up to 10 TeV. Above that energy, a 2-dimensional 
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FIG. 12: The estimated uncertainty in the atmospheric neutrino flux as a function of energy. Due 
to the large uncertainty in the prompt neutrino flux at greater than 10^ GeV, the total uncertainty 
rises sharply. 

fit of the lower energy values was again used as described above. The result was a lower 
atmospheric neutrino fiux prediction than the Barr et al. fiux. 

As described in Section IV At uncertainties in hadronic interactions and the cosmic ray 
and prompt neutrino fiuxes at high energy led to large total uncertainties in the atmospheric 
neutrino fiux. The estimated uncertainties are shown in Figure [121 

The atmospheric neutrino simulation was renormalized based on the experimental low 
energy data. The number of low energy conventional atmospheric neutrinos (second column) 
is added to the 4.0 prompt neutrinos predicted with the central prompt neutrino model. The 
total atmospheric background prediction before renormalization is shown in the third column 
of Table IXI[ Instead of renormalizing the simulation based on all events with A'ch < 100, 
the renormalization was only based on the region 50 < A^^ch < 100. Because of the difficulty 
of simulating events near the threshold of the detector, the atmospheric neutrino simulation 
did not faithfully reproduce the shape of the Ach distribution for the data at Ach below 50. 
Atmospheric neutrino models were scaled to match the 146 events seen in the experimental 
data for 50 < Ach < 100. The total number of high energy events predicted to survive 
the final energy requirement is shown before renormalization in the sixth column and after 
renormalization in the last column. 
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Barr et al. Max 
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253 
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14.5 


8.3 


Barr et al. 


194 
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0.74 


9.1 


10.3 


7.6 


Barr et al. Min 
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1.03 


4.9 


6.1 


6.3 


Honda et al. Max 


191 
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0.75 


9.3 


10.5 


7.9 


Honda et al. 


149 
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0.96 


6.4 


7.6 


7.3 


Honda et al. Min 


107 


111 


1.32 


3.4 


4.6 


6.1 



TABLE XI: Number of atmospheric neutrino events predicted by the Monte Carlo. Uncertainty in 
the high energy cosmic ray flux was incorporated into the maximum and minimum predictions. 
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Erratum: Multi-year search for a diffuse flux of muon neutrinos with AMANDA-II 

A search for TeV - PeV muon neutrinos with AMANDA-II data collected between 2000 
and 2003 estabhshed an upper hmit of E^$9o%c.l. < 7.4 x 10~^ GeV cm'^ s"^ sr"^ on the 
diffuse flux of extraterrestrial muon neutrinos with a $ oc spectrum between 16 TeV 
and 2.5 PeV. The upper limit calculation correctly included event simulations and remains 
as stated. However, the calculation of the detector's efficiency, which is based only on 
simulations, was incorrectly tabulated in an appendix and shown in a figure. The values 
were approximately a factor of ten too high, although the exact error varies in each bin. 
The correction has been apphed in the following tables and figure. The effective area is 
the equivalent area over which the detector would be 100% efficient for detecting neutrinos. 
The typical uncertainty on the effective area from simulation statistics is lowest between 10^ 
GeV and 10^ GeV (2%). The uncertainty increases to 6% at 10^ GeV and 5% around 10''' 
GeV. In the remainder of this document, the number of optical modules (OMs) triggered 
during an event is referred to as A^ch and cos(^t) refers to the cosine of the simulated (true) 
zenith angle of an event. The term angle- averaged indicates that results are averaged over 
9i between 100° and 180°. All other results reported in the paper, including the upper hmit, 
remain unchanged. 
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TABLE XII: Effective area as a function of the energy and zenith angle of the simulation for events 
in the final sample satisfying A^^^h ^ 100. 
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TABLE XIII: The angle-averaged neutrino effective area as a function of energy for events in the 
final sample satisfying N^ii > 100. 
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FIG. 13: Effective area for as a function of the true simulated energy at the Earth's surface 
in intervals of cosine of the true zenith angle, 9t- The angle-averaged effective area is represented 
by the solid black line. This calculation was based on the final event sample for events satisfying 
iVch > 100. 
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